Weak solution for the Hele-Shaw problem: viscous shocks and singularities 
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In Hele-Shaw flows a boundary of a viscous fluid develops unstable fingering patterns. At vanishing 
surface tension, fingers evolve to cusp-like singularities preventing a smooth flow. We show that the 
Hele-Shaw problem admits a weak solution where a singularity triggers viscous shocks. Shocks form 
a growing, branching tree of a line distribution of vorticity where pressure has a finite discontinuity. 
A condition that the flow remains curl-free at a macroscale uniquely determines the shock graph 
structure. We present a self-similar solution describing shocks emerging from a generic (2,3)-cusp 
singularity - an elementary branching event of a branching shock graph. 



1. Introduction Hele-Shaw flow describes a 2D vis- 
cous incompressible fluid with a free boundary at low 
Reynolds numbers. The fluid is either sucked out from 
a drain or driven to a drain by another, inviscid, incom- 
pressible liquid i]. Darcy law governs the viscous flow: 

j = -KVp, j = p v. (1) 

If density po and hydraulic conductivity K are constants, 
then pressure in the incompressible fluid p is a harmonic 
function, Ap = 0. At a drain (set at infinity), the fluid 
disappears with a constant flux Q = § j x d£. 

At vanishing surface tension (the case we consider), 
pressure is a constant along a boundary. Thus, in the 
fluid, pressure is a solution of the Dirichlet boundary 
problem. The boundary itself evolves in time. 

A compact form of the law involves only a boundary: 
normal velocity of a line element of the boundary d£ is 
proportional to its harmonic measure vd£ ~ dp (Har- 
monic measure is a distribution of a Brownian excur- 
sion (emanating from a drain) as it hits the boundary 
d/i = \Vp\d£). 

In this form, the Darcy law goes far beyond applica- 
tions to fluid dynamics 2] . It is closely related to a wide 
class of 2D growth and solidification processes such as 
DLA [3 , flows in granular media [4], visco-elastic flow 
[5], evaporative patterns of fluid monolayers [6], etc. 

Common patterns observed in these flows are charac- 
terized by intricate viscous fingering instabilities [7]. 

As such the problem is ill-defined: at a finite, criti- 
cal time fingers develop cusp singularities (points of in- 
finte curvature) [HHTU] , when the Darcy law stops making 
sense. Nevertheless, in many important physical prob- 
lems such as flows in granular media and solidification, 
which at large scales are described by Darcy law, flows 
are not limited and go over a singularity. Resolving finite- 
time singularities in a physical manner, and a description 
of a flow beyond singularities is a major long-standing 
problem in the field. 

An origin of singularities is the approximation of all 
physical parameters which could stabilize a flow at a mi- 
croscale being set to zero. While this is a valid approxi- 
mation at a smooth flow, at a singularity perturbations 



are singular, and an order of limits when different physi- 
cal parameters are brought to zero is essential. One per- 
turbation regularizing a flow is surface tension. However, 
interesting patterns are observed in processes where sur- 
face tension is either small, or does not exist at all, like in 
solidification and granular materials. This is the regime 
we are interested in. We assume that the surface tension 
(if any) must be set to zero first, before other parameters 
such as compressibility, are set to zero as well. 

In this Letter we present what we believe might be the 
solution to the problem of finite-time singularities in the 
Hele-Shaw problem. 

We impose the incompressibility and curl-free condi- 
tions at a large scale, but relax them at a microscale, 
setting surface tension to zero in the first place. Under 
this setting the Hele-Shaw problem admits a weak solu- 
tion: a singular cusp-like finger triggers shocks - lines of 
discontinuity for pressure. 

Shocks propagate through the fluid, forming a growing, 
branching tree. In this Letter we analyze a generic cusp- 
singularity giving rise of an elementary branching event 
of what will become a complicated degree-two tree of 
shocks. 

The solution describing local origin of viscous shocks 
or any further branching event of already existing shocks 
pattern is self-similar and does not depend on the details 
of the flow. 

A distinct feature of the Hele-Shaw flow is integrability 
[TTMl8j . Our weak solution is the only regularization of 
singularities which preserves the flow integrability. 

A comment is in order: the most common appear- 
ance of shocks in hydrodynamics are ultrasonic flows in 
compressible fluids at vanishing viscosity. Formally they 
are caused by inertial terms in hydrodynamic equations. 
Shocks in Hele-Shaw flow considered here are of a differ- 
ent nature. They occur at vanishing Reynolds numbers 
when inertial terms are negligible and viscosity is high. 
We refer to them as viscous shocks. 

2. Analytic and integrated forms of the Darcy Law We 
recall a description of a viscous finger in terms of a height 
function [16] . In Cartesian coordinates aligned with a fin- 
ger axis (FIG[l]), a fin ger is given by a graph y(x). Let 
X = x + iy be a complex coordinate of the fluid. The 
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FIG. 1: A growing finger in local Cartesian coordinates. The 
dashed line is the branch cut of the height function. 

height function Y(X, t) defined in the fluid is an analytic 
function whose boundary value on the fluid boundary is 
Y(X)x= x = y(x). Its discontinuity across the branch cut 
is discy|x=x = 2y(x). In terms of the height function 
the Darcy law reads: 

PqY = -d x <t>, (2) 

where the analytic function <j> = ip + ip is a complex 
potential of the flow, and ip is a stream function. 

An important physical characteristic is the (time- 
depedent) capacity of the flow. It is defined as C(t) = 
§ <t>d<f), where the integral goes around a drain. Ac- 
cording to Darcy law, the power required to drain the 
fluid N(t) =2¥k§ od P (jxd£) yields capacity N(t) = -C. 

Integrated form of Darcy law: The integral Q(X) = 
— i J poYdX gives yet another form of the Darcy law: 

ii = \4> = -p + iip, (3) 

Let us integrate ^ over a cycle B in the fluid: 

^Im j> dn = j> j x d£ = j> dip, (4) 
^-Re <£ dft = - I j • dt = - I dp = 0. (5) 

u-t J b J B J B 

The imaginary part measures a flux of fluid through the 
cycle, the real part measures circulation along the cycle. 
The latter vanishes. At infinity, i f df2 = Qt represents 
the mass of fluid drained up to time t. 

3. Singularities of the Hele-Shaw flow The pressure 
gradient is highest where the boundary curvature is large. 
Thus by (JlJ, growth velocity is largest (and increasing) 
at finger tips. It diverges at a critical time; the finger 
then becomes a cusp of type (2, 2/ + 1): y 2 ~ x 2l+1 jSJ- 

nni nsi lis] , fig.[i} 

Cusps of types (2, 4fc — 1) and (2, 4fc + 1) evolve differ- 
ently. In [17] it was shown that for the type (2,4k + 1), 
a new droplet of inviscid fluid emerges next to the fin- 
ger tip, before it evolves into a cusp. The fluid becomes 
multiply-connected, but evolution continues smoothly. 

No continuous evolution is possible for (2, 4k— 1) cusps, 
including the most generic cusp (2,3) - subject of this 
paper. 



4-. Hydrodynamics of a critical finger A finger ap- 
proaching the (2,3)-cusp is especially simple [15) . Fixing 
a scale and the origin, it is given by: 

-Y 2 = 4X 3 -g 2 X-g 3 = 4(X-e 3 )(X-e 2 ){X-e 1 ), (6) 

where gi and g 3 are real time dependent coefficients. One 
of the branch points, e 3 , may be chosen real. It is the tip 
of a finger. The other two are conjugated e\ — e 2 - The 
coefficient 172 = 4 ( | ei | 2 — e|) is determined by the drain 
rate, Q. We set the rate that g 2 — —I2t, and count time 
from a cusp event: flow is smooth at t < 0. 

In this normalization, capacity equals C = — <?2- Well 
before a critical time, e 3 < 0. Then at t < 0, g 2 > 
and Re e\ > 0: the two branch points are located in the 
fluid. This breaks the conditions of incompressibility, 
unless the roots coincide to a real double point: e\ = e 2 . 
Thus, g 3 = 8(-i) 3/2 . The height function at t < is a 
degenerate elliptic curve [15] : 

Y 2 = -4(X ~e{t))(X + e{t)/2) 2 , e(t) = -2a/=I (7) 

A finger becomes a cusp when the branch point e and 
double point — e/2 merge to a triple point. An important 
property of the critical finger Q is that it is self-similar: 

Y(X,t) = \t\ 3 ^Y{\t\- 1 / 2 X,l). (8) 

5. Weak form of the Hele-Shaw flow Once the flow 
reaches a cusp singularity, it is no longer governed by 
the differential form of the Darcy law. This situa- 
tion is typical for conservation laws of hyperbolic type 
dtu + d z f(u) = 0. There as well, smooth initial data 
develop into a shock at a finite time. Shocks occur to 
ill-defined conservation equations which arise as approx- 
imations of a well-defined problem. Adding a deforma- 
tion through terms with higher gradients prevents for- 
mation of singularities. However, a smooth solution of a 
deformed equation may become discontinuous as a de- 
formation is removed. It is then called a weak solu- 
tion [2^ • Validity of the differential equation on both 
sides of a shock determines a traveling velocity of a front 
V = jlscti " ^ ne Rankine-Hugoniot condition [2T] , 

Often, physical principles determine dynamics of 
shocks without specific knowledge of the deformation 
used, such that different deformations lead to the same 
weak solution. The best known example is the Maxwell 
rule determining the position of Burgers-type shock 
fronts EI]. 

Darcy law is a conservation law of hyperbolic type, 
and we adopt the same strategy. We will be looking for a 
weak solution of the Hele-Shaw problem when the Darcy 
law is applied everywhere in the fluid except on a moving, 
growing and branching graph T(t) of viscous shocks (or 
cracks), where pressure suffers a finite discontinuity. 

A few natural physical principles guide toward a 
unique weak solution. We give three equivalent formula- 
tions. 
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The first formulation is in terms of the height function, 
treated as a complex vector: 

- The canonically oriented (anti-clockwise) disconti- 
nuity of the height function's complex conjugate is 
parallel to the shock line directed such that: 



PodiscF|r = -2a£, a = p \Y\ > 



(9) 



where t is a unit vector canonically oriented along 
the shock line. 

An equivalent invariant formulation is in terms of the 
generating function (|3j|: 

- Discontinuity of £7 on shocks is imaginary; Re disc f2 
is increasing away from both sides of a shock, 



Re disc il I 



0, Ref7| x ^r>0. 



(10) 



The first condition in ( 10 ) is equivalently written as 
KeifYdX = for all cycles. Curves of this kind are 
called Krichever-Boutroux curves. They appeared in 
studies of Whitham averaging of non-linear waves [55] 
and asymptotes of orthogonal polynomials [23] . Neither 
appearance is coincidental. 

The third formulation is in hydrodynamics terms. 

6. Hydrodynamics of viscous shocks We require that 
the fluid is irrotational at macroscale, but allow vortic- 
ity at a microscale - i.e. at a scale set by a vanishing 
regularizing parameter. 

The first conditions of ( |9|T0 | insure that the fluid re- 
mains curl-free despite of shocks. Let us chose a contour 
B including a portion of a shock. Then ( 9|10 ) mean that 
the integral ^ still vanishes despite of discontinuities of 
fi and pressure. 

Let us study this condition in detail. The time 
derivative has two contributions: one from evolution 
of the height function Y(X,t) at a fixed X, another 
from motion of shocks. Denote the shock front veloc- 
ity by Vj_ (normal to the front, directed along the vec- 
tor n = — il). Then the time derivative in ^ reads 
Im [disc Y 1 + Vy (disc Y ■ Vj_)]r = 0, where Vy represents 
the derivative along the direction tangent to the front 
along the vector i. The first term is a jump of the veloc- 
ity of the fluid parallel to a shock —disc i?y . The second 
term is purely real. It equals Vy (crVjJ. Together, they 
yield to the condition: 

VyJ]_ + discjy =0, Jj_ = aV_L, j {] = p V {] . (H) 

The first term in this equation represents the transport 
of mass by a shock (normal to the shock), while the sec- 
ond is a circulation of the surrounding fluid flow. They 
compensate each other. 

This condition, derived solely from the requirement 
that a is real, suggests to interpret viscous shocks as 
a single layer of microscale vortices with a line density 
a and fixed orientation. Therefore, the weak solution 



describes a fluid with zero vorticity at a macroscale: vor- 
ticity concentrated in shocks is compensated by the cir- 
culation of fluid around shocks. 



Using Darcy law we replace the fluid velocity in (11) 
by — V and integrate (111 along the cut. We obtain 
the Rankine-Hugoniot condition: 



crVj_ = (disc p) n 



(12) 



The second set of conditions in ( 9|10 | [a > 0) imply 
that shocks move toward higher pressure, i.e, represent 
a deficit of fluid - cracks. This follows from the curl-free 
condition and is consistent with an assumption that at 
infinity (at a remote part of a finger, and at a drain), the 
flow is not affected by shocks. 

Consider the integral around the drain i f dQ, = Qt. 
Before the transition, a contour of integration can be 
smoothly deformed to the fluid boundary. It yields the 
total mass 2p § y(x)dx = po JJ dydx. After the tran- 
sition, the integral acquires a contribution from shocks: 
— <f> r ad£. Viscous shocks can be seen as a single layer 
mass deficit: density of the fluid is no longer constant, 
p{X) = po — 5r(X)a(X) (here 5r is the delta- function on 
shocks). 



The Rankine-Hugoniot condition (12), curl- free condi- 
tion ^ and the differential form of the Darcy equation 
([2]) combined give the weak form of the Darcy law. Re- 
markably, these conditions uniquely determine a shock 
pattern. Shock graphs with one and two branching gen- 
erations are represented in FIGs. 2,3. We will show how 
they work by an analysis of the (2,3)-cusp . 

7. Self-similar weak solution: beyond the (2,3)-cusp 
In the remaining part of the Letter we briefly describe a 
solution to the most generic singularity, (2,3). 

Before becoming cusps, fingers are described by Q. 
The height function is a self-similar degenerate elliptic 
curve. Two branch points located in the fluid coincide 
to a double point. At a critical time, the tip meets the 
double point. After the critical time, we look for a weak 
solution, allowing pressure to have finite discontinuities 
on some curves. We assume that infinity is not affected 
by the transition, such that the height function is still an 
elliptic curve (JsJ) with g 2 = — 12t, but now g 2 < 0. Wc 
have only to find g 3 from the condition (10). Solution 



is again self-similar, g 3 (t) = g 3 (l)t 3 / 2 . Once g 3 (l) ^ 8, 
the curve is not degenerate anymore. The double point 
splits into two branch points, e\ 7^ e 2 . They appear in 
the fluid as endpoints of shocks. 

Since branch points are simple, all quantities have op- 
posite signes on opposite sides of shocks: disc SI = 2f2|r- 



Then condition ( 10 1 means that Re Q vanishes on shocks. 
Since the endpoints belong to the fluid, by virtue of (|3j, 
pressure also vanishes. The branch point e 3 is a tip of 
a finger where pressure obviously vanishes. This is the 
governing condition: 

p(ei, a ) =W(ei, 2 ) = 0. (13) 
The scaling property (Tsl) is sufficient to express the po- 
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FIG. 2: The equi-pressure lines and flow lines (dashed) before 
the transition (left), and after the transition (right). Pressure 
gets larger as the shade gets darker. Before the transition 
62 — ei is the double point. The thick lines (on the right 
panel) connecting the branch points are the shocks. The or- 
ange and green dashed lines are not admissible the level lines 
of Re n. The fluid flows to the lighter region toward low pres- 
sure. Shocks moves to toward the darker region (higher pres- 
sure). A bright dotted line is the zero-pressure line. The ar- 
rows are moving directions of the shock and the branch point. 
On the left of the zero-pressure line the finger expands, on the 
right the finger retreats. 



tential cj) and Q through elliptic integrals: 

' ' Je 3 \ 2g 2 jY(X)- 

r x oi 
fl(X) = -i YdX = (XY-2t(j)). (15) 

Je 3 5 



(14) 



The governing equation ( 13 1 is then expressed through 
complete elliptic integrals K and E [23] : 

(16m 2 - 16m + l)E{m) = (8m 2 - 9m + l)K(m) , (16) 
where m = \ + 1^^= and f ^ = z Jj7 3 %t\ - 



Solution of (16) is transcendental [24] 
ei,2 1 



± ih. 



e3 



3.24638225374 . 



(17) 



It determines the moving ends of shocks. 

8. Discontinuous change of capacity The genus tran- 
sition is summarized by an abrupt change of g?,\t=±i- We 
express this fact in normalization-independent terms us- 
ing time derivative of the capacity. The ratio 



Ct>0 _ ^53.-1/2 

Ct<o ~ 2 92 



lim 



0.91522030388 . 



(18) 



is a unique universal number describing the transition. 

9. Flow and shocks. Shocks are anti-Stokes lines of 
f2, i.e., zero- level lines of Refi selected by the admissi- 
bility condition ([9j [l0]) fi(X)|x->-r > 0. Refi vanishes 
at Im (XY) = 2tp(X). There are a total of seven anti- 
Stokes lines connected at branch points. They are tran- 
scendental and computed numerically, FIG. [2] 




FIG. 3: A numerically computed graph of viscous shocks with 
two generations of branching. 



Among the seven anti-Stokes lines, only three lines 
T3,r2 and Ti obey the second condition (10). Line Ti is 
the boundary of the finger. Lines r3,r2 are shocks. 

Selection works as follows. A remote part of the finger 
(argX = 7T, \X\ — > oo) is not affected by the transition. 
This selects a branch of fl ~ |X 5 / 2 such that on the 
upper side of T%, Kefl > 0. and therefore in both sectors 
{IV, III} where |argA" - tt| < 2?r/5 at a large X. Then 
signs of Re are opposite on both sides of T4 , T 5 , r 6 , , , 
and signs are positive on both sides of Fi,r2,r3, as re- 
quired by ( 10 ). 



The Rankinc-Hugoniot conditions (12 and |9]) give the 
velocity of shocks. Noticing that Im XY = crX±, where 
X± is a projection of a vector-coordinate of a shock to a 
direction normal to the shock, we get Vj_ — 

Already scaling yields that shocks push the fluid away 
faster (^t 5 ^ 4 ) than it is absorbed by the drain ~t. There- 
fore, the finger retreats v(e$) = £3 < 0, smoothing the 
tip. Indeed, at the endpoints ^(ei^) is purely imagi- 
nary. Then by virtue of |9| Im f2(ei,2) = J^' 2 ad£ is a 
mass deficit accumulated on shocks. Eqs.( |15|17 ) yield: 



ff<W=-|V(ei, 2 )| 
5 



- (6.34513... )t 



5/4 



(19) 



Since the boundary of the fluid moves towards lower pres- 
sure, pressure is positive close to sides of the tip, but 
remains negative around a distant part of the finger. 
Therefore, in the fluid there are lines of zero pressure 
(the bright dotted line in FIG. [2]). They emanate from 
the end points crossing Ti normally at a point x = —rjt^' 2 , 
where velocity vanishes, FIG. [2] 

Finally, we list angles of zero-pressure line, the shock 
line and velocity at e\. Relative to the real axis, they are 
respectively 0.235061tt, 0.3792582tt, and 0.451357tt. At 
e3 the angle between shocks is 27r/3. 

10. Branching tree Evolution through a (2,3)-cusp 
gives rise to a shock tree. It grows and keeps branching 
further. An interesting branching tree emerges, FIG. [3] 
We do not know its global structure, but we do know 
that every generic branching event is locally identical to 
the transition we just described. It will be interesting 
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to study whether a developed shock's graph with a large 
number of branches exhibits a scale invariance. 

11. Viscous shocks may have different meanings de- 
pending on experimental settings. In a Helc-Shaw cell, 
viscous shocks are narrow channels where an inviscid liq- 
uid is compressed and sheared. In a visco-ellastic fluid 
and in granular materials viscous shocks are cracks, etc. 
It will be very interesting to see an experimental realiza- 
tion of the branching tree of shocks in viscous flows at 
small Reynolds numbers. 
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